Dynamical configuration interaction: Quantum embedding that combines wave functions and Green's functions

We present the concept, derivation, and implementation of dynamical configuration interaction, a quantum embedding theory that combines Green's function methodology with the many-body wave function. In a strongly-correlated active space, we use full configuration interaction (CI) to describe static correlation exactly. We add energy dependent corrections to the CI Hamiltonian which, in principle, include all remaining correlation derived from the bath space surrounding the active space. Next, we replace the exact Hamiltonian in the bath with one of excitations defined over a correlated ground state. This transformation is naturally suited to the methodology of many-body Green's functions. In this space, we use a modified $GW$/Bethe-Salpeter equation procedure to calculate excitation energies. Combined with an estimate of the ground state energy in the bath, we can efficiently compute the energy dependent corrections, which correlate the full set of orbitals, for very low computational cost. We present dimer dissociation curves for H$_2$ and N$_2$ in good agreement with exact results. Additionally, excited states of N$_2$ and C$_2$ are in excellent agreement with benchmark theory and experiment. By combining the strengths of two disciplines, we achieve a balanced description of static and dynamic correlation in a fully ab-initio, systematically improvable framework

A quantum embedding theory in the screened Coulomb interaction: Combining configuration interaction with GW/BSE

We present a new quantum embedding theory called dynamical configuration interaction (DCI) that combines wave function and Green's function theories. DCI captures static correlation in a correlated subspace with configuration interaction and couples to high-energy, dynamic correlation outside the subspace with many-body perturbation theory based on Green's functions. In the correlated subspace, we use a wave function description to avoid embedding the two-particle vertex, which greatly simplifies the frequency structure of the embedding. DCI takes the strengths of both theories to balance static and dynamic correlation in a single, fully ab-initio embedding concept. We show that treating high-energy correlation up to the $GW$ and Bethe-Salpeter equation level is sufficient even for challenging multi-reference problems. Our theory treats ground and excited states on equal footing, and we compute the dissociation curve of N$_2$, vertical excitation energies of N$_2$ and C$_2$, and the ionization spectrum of benzene in excellent agreement with high level quantum chemistry methods and experiment

Piecewise linearity in the GWGW approximation for accurate quasiparticle energy predictions

We identify the deviation from the straight line error (DSLE) -- i.e., the spurious non-linearity of the total energy as a function of fractional particle number -- as the main source for the discrepancy between experimental vertical ionization energies and theoretical quasiparticle energies, as obtained from the $GW$ and $GW$+SOSEX approximations to many-body perturbation theory (MBPT). For self-consistent calculations, we show that $GW$ suffers from a small DSLE. Conversely, for perturbative $G_0W_0$ and $G_0W_0$+SOSEX calculations the DSLE depends on the starting point. We exploit this starting-point dependence to reduce (or completely eliminate) the DSLE. We find that the agreement with experiment increases as the DSLE reduces. DSLE-minimized schemes, thus, emerge as promising avenues for future developments in MBPT

Thermodynamic equilibrium conditions of graphene films on SiC

First-principles surface phase diagrams reveal that epitaxial monolayer graphene films on the Si side of 3C-SiC(111) can exist as thermodynamically stable phases in a narrow range of experimentally controllable conditions, defining a path to the highest-quality graphene films. Our calculations are based on a van der Waals corrected density functional. The full, experimentally observed (6 sqrt(3)x 6 sqrt(3))-R30 supercells for zero- to trilayer graphene are essential to describe the correct interface geometries and the relative stability of surface phases and possible defects

Exact-exchange based quasiparticle energy calculations for the band gap, effective masses and deformation potentials of ScN

The band gaps, longitudinal and transverse effective masses, and deformation potentials of ScN in the rock-salt structure have been calculated employing G0W0-quasiparticle calculations using exact-exchange Kohn-Sham density functional theory one-particle wavefunctions and energies as input. Our quasiparticle gaps support recent experimental observations that ScN has a much lower indirect band gap than previously thought. The results are analyzed in terms of the influence of different approximations for exchange and correlation taken in the computational approach on the electronic structure of ScN.Comment: Comments: 9 pages including 5 figures; related publications can be found at http://www.fhi-berlin.mpg.de/th/th.htm

Self-Consistent Green Function Embedding for Advanced Electronic Structure Methods Based on a Dynamical Mean-Field Concept

We present an embedding scheme for periodic systems that facilitates the treatment of the physically important part (here the unit cell) with advanced electronic-structure methods, that are computationally too expensive for periodic systems. The rest of the periodic system is treated with computationally less demanding approaches, e.g., Kohn-Sham density-functional theory, in a self- consistent manner. Our scheme is based on the concept of dynamical mean-field theory (DMFT) formulated in terms of Green functions. In contrast to the original DMFT formulation for correlated model Hamiltonians, we here consider the unit cell as local embedded cluster in a first-principles way, that includes all electronic degrees of freedom. Our real-space dynamical mean-field embedding (RDMFE) scheme features two nested Dyson equations, one for the embedded cluster and another for the periodic surrounding. The total energy is computed from the resulting Green functions. The performance of our scheme is demonstrated by treating the embedded region with hybrid functionals and many-body perturbation theory in the GW approach for simple bulk systems. The total energy and the density of states converge rapidly with respect to the computational parameters and approach their bulk limit with increasing cluster (i.e., unit cell) size

Renormalized Second-order Perturbation Theory for The Electron Correlation Energy: Concept, Implementation, and Benchmarks

We present a renormalized second-order perturbation theory (rPT2), based on a Kohn-Sham (KS) reference state, for the electron correlation energy that includes the random-phase approximation (RPA), second-order screened exchange (SOSEX), and renormalized single excitations (rSE). These three terms all involve a summation of certain types of diagrams to infinite order, and can be viewed as "renormalization" of the 2nd-order direct, exchange, and single excitation (SE) terms of Rayleigh-Schr\"odinger perturbation theory based on an KS reference. In this work we establish the concept of rPT2 and present the numerical details of our SOSEX and rSE implementations. A preliminary version of rPT2, in which the renormalized SE (rSE) contribution was treated approximately, has already been benchmarked for molecular atomization energies and chemical reaction barrier heights and shows a well balanced performance [Paier et al, New J. Phys. 14, 043002 (2012)]. In this work, we present a refined version of rPT2, in which we evaluate the rSE series of diagrams rigorously. We then extend the benchmark studies to non-covalent interactions, including the rare-gas dimers, and the S22 and S66 test sets. Despite some remaining shortcomings, we conclude that rPT2 gives an overall satisfactory performance across different chemical environments, and is a promising step towards a generally applicable electronic structure approach.Comment: 16 pages, 11 figure

Self-interaction in Green's-function theory of the hydrogen atom

Atomic hydrogen provides a unique test case for computational electronic structure methods, since its electronic excitation energies are known analytically. With only one electron, hydrogen contains no electronic correlation and is therefore particularly susceptible to spurious self-interaction errors introduced by certain computational methods. In this paper we focus on many-body perturbation-theory (MBPT) in Hedin's GW approximation. While the Hartree-Fock and the exact MBPT self-energy are free of self-interaction, the correlation part of the GW self-energy does not have this property. Here we use atomic hydrogen as a benchmark system for GW and show that the self-interaction part of the GW self-energy, while non-zero, is small. The effect of calculating the GW self-energy from exact wavefunctions and eigenvalues, as distinct from those from the local-density approximation, is also illuminating

Density-functional Theory for f electron Systems: the {\alpha}-{\gamma} Phase Transition in Cerium

The isostructural {\alpha}-{\gamma} phase transition in cerium is analyzed using density-functional theory with different exchange-correlation functionals, in particular the PBE0 hybrid functional and the exact- exchange plus correlation in the random-phase approximation [(EX+cRPA)@PBE0] approach. We show that the Hartree-Fock exchange part of the hybrid functional actuates two distinct solutions at zero temperature that can be associated with the {\alpha} and {\gamma} phases of cerium. However, despite the relatively good structural and magnetic properties, PBE0 predicts the {\gamma} phase to be the stable phase at ambient pressure and zero temperature, in contradiction with low temperature experiments. EX+cRPA reverses the energetic ordering, which emphasizes the importance of correlation for rare- earth systems